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Abstract 

We consider proton stability based on inspired extra U(l) model with 5*4 flavor symmetry. In 
this model, a long life time of proton is realized by the flavor symmetry. One of the interesting effects 
of flavor symmetry is that the proton decay widths of p fi + X are suppressed by cancellation. This 
^ C"| suppression mechanism is important in the case that Yukawa coupling constants are hierarchical. Our 

O j model predicts p — > e + K° has larger decay width than that of p — > pL + K°. 
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1 Introduction 



Supersymmetry is an elegant solution of hierarchy problem of Standard Model (SM) [1], however, a simple 
supersymmetric extension of SM suffers from non-conservation of baryon number and /it-problem. Therefore 
we must introduce new symmetry such as R-parity in MSSM, to suppress proton decay operators and fj,- 
term. One of the solutions of /i-problcm is given by introducing extra U(l) gauge symmetry [2]. In this 
frame work, several new superfields such as singlet S, exotic quarks G,G C , must be introduced to cancel 
gauge anomaly. This is the elegant solution of /i-problem, however proton instability is not solved because 
the baryon number violating interactions in superpotential of MSSM are replaced by single exotic quark 
interactions. 

In the superpotential, there is no obvious distinction between baryon number violating trilinear terms 
and Yukawa interactions. Therefore it is natural to introduce flavor symmetry in suppressing baryon number 
violating operators. Especially non-Abelian discrete symmetry is good candidate for the flavor symmetry, 
because large mixing angles of Maki-Nakagawa-Sakata (MNS) mixing matrix may be explained simultane- 
ously, and more simply, non-Abelian symmetry can be the reason why generation exists. At previous work, 
we explained S4 x Z2 flavor symmetry not only realizes maximal mixing angle 623 but also suppresses pro- 
ton decay based on SU(3) C x SU(2)w x U(l)y x U(l)x x U(l)z gauge symmetry [3]. As the suppression 
mechanism of proton decay is complicated and model dependent, we give more detailed estimation of proton 
life time and investigate several flavon sectors in this paper. 

In the several models of the non-Abelian discrete symmetries, although the direction of vacuum expec- 
tation values (VEVs) of flavons is crucial, the direction is selected by hand adding the explicit symmetry 
breaking terms. As such procedure is unnatural, we start from exactly flavor symmetric theory in the es- 
timation of proton life time. The spontaneous flavor symmetry breaking realizes special VEV direction of 
flavons, which affects proton life time significantly. If the Yukawa hierarchy is realized by Froggatt-Nielsen 
mechanism, as p — > e + X are suppressed due to the small coupling constants, p — > fi + X may dominate proton 
decay width. In our model, as p — > [i + X are suppressed by cancellation, p — > e + K° dominates the proton 
decay width. 

This paper is organized as follows. In section 2, we estimate proton life time based on gauge non-singlet 
flavon model. In section 3, we modify the flavon sector by adding gauge-singlet flavon and Froggatt-Nielsen 
flavon. In section 4, we eliminate gauge non-singlet flavon and construct Dirac neutrino model. Finally we 
give conclusion of our analysis in section 5. 

2 S4 x Z2 flavor symmetric extra U(l) model 

At first we explain the basic structure of our model. We extend the gauge symmetry to G2 = SU(3) C x 
SU(2)w x U(1)y x U(l)x x U(l)z which is the maximal subgroup of Eq, In order to cancel gauge 
anomaly, we must add new superfields, such as SM singlet S, exotic quark G, G c (hereafter we call them 
g-quark) and right handed neutrino (RHN) N c . We can embed these superfields with MSSM superfields 
Q, U c , D c , L, E c , H u , H D into 27 of E 6 [4]. As the singlet S develops VEV and breaks U(l) x gauge sym- 
metry, O (TeV) scale /i-term is induced naturally. In order to break U(l)z and generate a large Majorana 
mass of RHN, we add SM singlet $, $ c . The gauge representations of superfields are given in Table 1 [3]. 
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Table 1: G 2 assignment of fields. Where the x, y and z are charges of U(l)x, U(l)y and U(l)z, and Y is 
hypcrcharge. R-parity R = cxp [§jj(3a; — 8y + 15z)] is unbroken. 

Under the gauge symmetry given in Table 1, the renormalizable superpotential is given by 

W = Y U H U QU C + Y D QD C H D +Y E H D LE C + Y N H U LN C + Y M <S>N C N C + XSH U H D + kSGG c 
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+ Af$$ c + Y QQG GQQ + y udg G c U c D c + Y EUG GE C U C + Y QLG G C LQ + Y NDG GN C D C . (1) 

In this superpotential, unwanted terms are included in the second line. The first term of the second line is 
the mass term of singlets <5, $ c which prevent singlets from developing VEVs. The other five terms of the 
second line are single g-quark interactions, which break baryon and lepton number and induce rapid proton 
decay. In the first line, we must take care of the flavor changing neutral currents (FCNCs) induced by extra 
Higgs bosons [5]. Therefore the superpotential Eq.(l) is not consistent at the present stage. 

In order to stabilize proton, we introduce S4 x Z2 flavor symmetry. If we assign G, G c to S4 triplet and 
quarks and leptons to doublet or singlet, the single g-quark interaction is forbidden. However, as the g-quark 
must never be stable from phenomenological reason, we assign $ c to S4 triplet to break the flavor symmetry 
slightly. In this case, as 4>, $ c play the role of flavons, we call them gauge non-singlet flavons. In order to 
realize the large mixing angle of #23 in the MNS matrix and suppress the Higgs-mcdiatcd FCNCs, we assign 
the superfields in our model as given in Table 2 [6]. 

In the non-renormalizable part of superpotential, the single g-quark interactions which contribute to the 
g-quark decay are given as follows 
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Table 2: S4X Z2 assignment of superfields (Where the index i of the S4 doublets runs i = 1,2, and the index 
a of the S4 triplets runs a = 1, 2, 3.) 



2.1 Higgs sector and hidden sector 

Under the flavor symmetry given in Table 2, the superpotential of Higgs sector is given by, 

W H = \iS 3 (H?H? + H^Hi > ) + X 3 S 3 H^Hi ) 

+ XiH^SiH? + S 2 H 2 D ) + \ 5 (S 1 H¥ + S 2 h¥)H?, (3) 

where one can take, without any loss of the generalities, Ai 3 4 5 as real, by redefining four arbitrary fields of 
{Si, 5*3, Hf , ,H[ > , H®}. As only ^4 singlets , and S3 couple to quarks and g-quarks respectively, 
they behave like MSSM Higgs and SM singlet respectively. Through the renormalization group equations, 
the squared masses of H^,S3 become negative, they develop VEVs and break gauge symmetry. As the 
result, the A-term A3S3H3 enforces developing VEV. However, S4 doublets do not develop VEVs. 
To generate the VEVs of S4 doublets, we must add flavor breaking squared mass terms. The origin of flavor 
breaking terms is discussed below. Note that there is accidental 0(2) symmetry induced by the common 
rotation of the S4 doublets in Wh ■ 

As it is thought that the scalar squared masses are induced as the result of SUSY breaking in hidden 
sector, we assume flavor symmetry is broken at the same time. We assume hidden sector is described by 
flavor symmetric extension of O'Raifeartaigh model [7]. We introduce gauge singlet A, B + , Bi,C+,Ci and 
assign Z'^ charges to them to separate hidden sector from observable sector. We assume U(1)r symmetry is 
hold at the limit of infinite Planck scale, Mp — > 00. The representations of hidden sector superfields under 
the S4 x Z2 x Z' 2 x U(1)r symmetry are given in table 3. 

Under the symmetry given in Table 3, the superpotential of hidden sector is given by, 

^hidden = -M 2 A + m + B + C + + m(B 1 C 1 +B2C2) + ^X + ACl + ^XA(C 2 1 +Ci). (4) 
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Table 3: S4 x Z 2 ( 4 ) x Z 2 x [/(l)i?. assignment of superficial (Where the index i of the £4 doublets runs 
1 = 1,2.) 
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SUSY is spontaneously broken since it is impossible to satisfy the equations Fa = and F B+ = -Fs, = at 
the same time. Notice that the flavor symmetry 64 x Z 2 is also broken spontaneously. As the superpotential 
has accidental 0(2) symmetry, the direction of £4 doublet F-term defined by, 



F Bl = Fc B , F B2 = F Sl 



(12) 



is described by free parameter 9 B - Although the spontaneous breaking of 0(2) results the appearance of 
Nambu-Goldstone boson (NGB), we assume this NGB does not cause any problem because the interactions 
of the NGB with observable sector particles are suppressed by Planck scale. 

The SUSY breaking in hidden sector is mediated to observable sector by gravity through the non- 
renormalizable terms in Kahler potential as given by, 



J ff 2 S 2 )( J ffiB 1 + H 2 B 2 y + (b H B+H 3 {H 1 B 1 + H 2 B 2 y + /i.e.)], (13) 

L p 

where H = H u , H D , S, and flavor symmetric terms and the contributions to other supcrficlds arc omitted. 
These terms induce scalar squared mass terms as follows 



V F b = mgjjAcBHx + s B H 2 \ 2 + [m 2 BH2 H 3 (c B Hi + s B H 2 y + h.c 



If we substitute the VEV (Hi) = vch, (H 2 ) = vsh, (H3) = v' for H a , then we get 
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At second line, we simplified the equation, because we are interested only in direction 9h- As we can change 
the sign of b by the redefinition of the sign of v, we can define b > without loss of generality. The minimum 
of potential Vfb is classified as follows 



a > 0, b > 1 
a > 0, b < 1 
a < 



?H = 



arccosfr, 

7T, 



(16) 
(17) 
(18) 



from which one can see that the angle 9h is controlled by free parameters a and b for the given Ob- In this 
section, we assume the condition Eq.(16) is satisfied for H u , H D , S and select the common VEV direction 
9 B . In this direction, the flavor symmetric part of Higgs potential has accidental 0(2) symmetry and do not 
depend on Oh and 9h that are fixed by Vfb- Therefore we can not assume Vfb as perturbation, even if 
flavor breaking parameter is small. 
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2.2 Flavon sector 

The superpotential of flavon sector is given by. 
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2M P 

and the flavor symmetric part of potential is given by, 

V = -mfl^l'+miO^ + l^n+mlO^ + l^^ + l^l 2 ] 

f $ 2 [(^) 2 + (*2) 2 + (1>3) 2 ] 
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r ^$ 3 {V3$ i [(*§) 2 - (*I) 2 ] + $ 2 [(<^) 2 + ($^) 2 - 2($J) 2 ] } 
+ F-term + D-tcrm. (20) 

As this potential does not have accidental continuum symmetry, flavor breaking terms can be assumed as 
perturbation. We assume negative squared mass of $3 pulls the trigger of U(l)z breaking and $3 and 
<I>q develop VEVs along the D-flat direction. For the VEV of $j, there are two possibilities in the flavor 
symmetry breaking. The one of them is 63 symmetric vacuum defined as, 

$3 = ^, $, = 0, ^ = ^(1,1,1), (21) 

and the other is S3 breaking vacuum defined as, 

$ 3 = Vc, $i=V8(0,i), *| = ^(0,1,1). (22) 

Another degenerated vacuums are given by acting $4 translations on these VEVs, respectively. Note that 
spontaneous breaking of discrete symmetry causes domain wall problem. We assume flavor symmetry is not 
recovered in reheating era after inflation [8]. As the vacuum defined in Eq.(22) is not phcnomenologically 
allowed because g-quark gi, g\ become stable, we select the vacuum defined in Eq.(21) by tuning m| large 
enough. 

From the simplified potential 

^-mWl*f + lJl*T. (23) 



M P 



the size of flavon VEV is estimated as follows 



(24) 

If we put y* ~ 0.01 and msuSY ~ 10 TeV, then we get V ~ 10 12 GeV that is a favorable value to satisfy 
the constraint for g-quark and proton life time at the same time [9] . 
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2.3 Quark and Lepton sector 

As the flavor symmetry reduces the number of free parameters drastically, we can decide the value of free 
parameters by very few assumptions. For the quark sector, superpotential is given by 



Wq = Y» b HgQ a U§ + Y»H»Q a D 



(25) 



As the extra Higgs H^Hf do not couple to quarks, Higgs mediated flavor changing neutral currents are 
not induced. 

For the lepton sector, superpotential is given by 
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Without any loss of generalities, by the field redefinition, we can define Y E 2 3 ,Yi 3 ,Y 2 N 4 are real and non 
negative and Y 2 M , Y 3 N are complex. We tune the angle 9b = & 2 : 
follows 
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and define the mass parameters as follows [10] 

Mj = YfV, M 3 
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m\ = Yfv d , m l 2 ■ 

With these parameters, the mass matrices are given by 
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Due to the seesaw mechanism, the neutrino mass matrix is given by 
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The charged lepton mass matrix is diagonalized as follows 

V}MtMtVi = diag(m 2 e ,ml,m 2 T ) = ((m 2 ) 2 , (m 3 ) 2 , (m[) 2 ), 
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and those of the light neutrinos are given by 

V^MyVv = diag^^W^e^+^m^m^), 

Vu 

From the above equations, the MNS matrix is given by 

1 
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As the value of 6b is correctly tuned, experimental value of 023 is realized. From the experimental constraints 
fill. 
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If we put the size of VEVs of Higgs fields as follows 

Vu = 10, v' u = 155.3, v d = 2.0, v' d = 77.8 (GcV), 
we can fix Yukawa coupling constants as follows 

Yf = 0.875, Y 3 E = 5.15 x 10" 2 , Y 2 E = 6.25 x 10" 6 , 
where we use the charged lepton masses [12] as follows 

m[ = 1.75 GcV, m\ = 487 kcV, m\ = 103 McV. 
For the neutrinos, if we put 

V = 10 12 GcV, Yj M = Y 3 M = 1, 
and assume = 0, all Yukawa coupling constants are determined as follows 
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pi = 1.795 x 10~ 2 cV, p\ = 15.51 x 10~ 2 eV, p\ = 13.79 x 10~ 2 eV, 
m Vl = 5.240 x 10~ 2 eV, m U2 = 5.312 x 10~ 2 eV, m V3 = 1.795 x 10~ 2 eV, 
m y 2 = 4.24 GcV, m\ = 12.45 GcV, m\ = 11.74 GcV, 
Y 9 N = 0.424, Yo N = 0.080, Yf = 1.17. 



(44) 



For the lepton sector, there is no flavor changing process as same as quark sector discussed above. Considering 
the interactions 
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(45) 



/i — > e7 and r — > processes are forbidden because e does not couple to /z and r. Furthermore, due to 
the unbroken 5*2 symmetry such as l fl —> — / /J ,/i c — » — p, c , (HP,HP) —> (HP , iff), r — >• ^7 process is also 
forbidden. 
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2.4 g-quark sector 

As the masses of g-quarks are given by 

W G = kS^GxGl + G 2 G C 2 + G 3 G C 3 ), (46) 

the g-quark mass matrix is proportional to unit matrix. For the scalar g-quarks, as the contribution from 
soft flavor breaking term should be added, the degeneracy of masses may be broken. However, such effects 
can be assumed as perturbation. It is thought that dark matter does not have strong interaction, g-quark 
should not be stable if the reheating temperature is higher than g-quark mass. Under the symmetry defined 
in Table 1 and Table 2, g-quarks can decay through the non-renormalizable terms as follows 



W B = 
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-^Yg DG ^[G^\ + G&l + Gt^U^Dt 
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+ ^r a LdQG $3[^£i(G^^-G§$§) + L 2 (G^- + G§^-2G 1 $^)]Q a . (47) 

Here we estimate the g-quark life time by the Y EUG interaction. If we put 

Y EUG = Y 2 EUG = Y EUG = 1, (48) 

then the interaction is given by 

Cg = {A %mT El9x{ul + u% + u§) ' (49) 

where E 2 is right handed selectron. The renormalization factor (A E ,p G )s is evaluated by the RGE 

„ .dlnY EUG 16 . 

4t — jt = — 5-«a, (50) 

a In fi 6 

and given by 

/ „ , x 4a,/37r /n .„ in i8\ 0.05008 

{ASO)s _ ^ .(^-) -«. T , (.., 

where only QCD correction is accounted. This approximation is not bad because the beta function of the 
coupling constant of strong interaction g s vanishes at 1-loop level in our model, which makes the contribution 
of a s dominant in the RGE of Y EUG . 

Using the interaction Eq.(49), we calculate the g-quark decay width. For simplicity, we assume u c a in 
Eq.(49) are mass eigenstates. Requiring the life time of g-quark is shorter than 0.1 sec (otherwise the success 
of BBN is spoiled [13]) as follows 

rC9l) -H 3MJ, ) 16^ > 0T^' (52) 

we get 

_i&_(0 >2 , 25xlO - (53) 
which bounds the VEV size of flavon from below. 
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Finally, using the interaction Y®Q G — Y EUG , we estimate the proton decay width. Integrating out the 
scalar g-quarks, we get the effective four-Fermi interactions as follows 
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A RF = 159.5. 

In the quark mass basis, Eq.(54) is rewritten as follows 
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For simplicity, we put 
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then the proton decay width is given by 
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If we put 



^ = 0.47, £> = 0.80, a p = -0.012 GeV 3 , /„ = 130 MeV, [17] 
m n o = 135 MeV, m p = 940 MeV, [11] 



then we get 
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From the experimental bound r(p — > ir° + e + ) > 1600 x 10 30 [years] [11], the VEV size of flavon is bounded 
from above as follows 



V \ 4 /1000 GeV x 



M P ) V M c 



< 2.60 x 10~ 54 . 



(66) 
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Hereafter we assume the approximation M g = Mq is held for simplicity. From Eq.(53) and Eq.(66), the 
allowed region for V is given by (sec Fig.l) 

2.25 x 10- r i000GeV ) < f-LV < 1.61 x 10-" ( MG V. (67) 
V Ma J \MpJ VIOOOGeVy v ' 

This inequality holds when the mass bound, 

Mq > 2.41 TcV, (68) 

is satisfied. 
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Figure 1: Mq versus V: The pink region comes from the constraint of the life time of the g-quark, which 
should be less than 0.1 sec. The green region comes from the constraint of the proton stability. The black 
region is allowed by the both constraints. The heavier of Mq, the wider the allowed region is. 



Before ending this section, we discuss the unsatisfactory point of this model. Considering the mass 
spectrum of the quarks and charged leptons, it is expected that the trilinear coupling of first generation is 
multi suppressed by the suppression mechanism of Yukawa couplings and 54 symmetry. If it is true, as the 
proton decay width accommodates another suppression factor, the experimental verification of proton decay 
seems to be impossible. Therefore the suppression by the gauge non-singlet flavon may be too strong. To 
improve this point, we modify the flavon sector in next section. 



3 S4 x x Zg flavor symmetric model 

In this section we introduce Froggatt- Nielsen mechanism to explain Yukawa hierarchy [18], and the flavon 
sector is modified as follows. To realize O(10 2 ) hierarchy, we introduce Z$ symmetry and add gauge and 
S4 singlet X as Froggatt-Nielscn (FN) flavon. To weaken the over suppression of trilinear terms, we replace 
the flavon by S4 singlet $ c and gauge singlet T which is assigned to S4 triplet. To forbid renormalizablc 
terms of T, the Z2 symmetry is replaced by Z4. The flavor representations of supcrficlds arc given in Table 
4. 



3.1 Flavon sector 



The leading terms of flavons are given as follows 
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W 7 



W T - 
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(70) 
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Table 4: S4X Z4X Zg assignment of superfields (Where the index i of the S4 doublets runs i = 1,2, and the 
index a of the S4 triplets runs a = 1,2,3.) 



6M P 



X 9 , 



2Mi 



-y*($$ c ) 2 . 



(71) 
(72) 



The VEV size of gauge non-singlet is estimated by Eq.(24). Now we change the value of V — ( < f>) to 10 11 
GeV which is given by naive estimation when we put Y* ~ 1. This affects neutrino Yukawa couplings given 
in Eq.(44) as follows 



1 2 



1 2 



0.134, Y, 



N 



*3 



TO V10 
The VEV size of FN flavon is estimated as follows 



0.025, V 



N 



V N 



0.370. 



U) 



/ msusY 
\ M P 



10" 



(73) 



(74) 



S4 symmetric part of potential of T is given by 

V T = -m 2 (|T!| 2 + |T 2 | 2 + |T 3 | 2 ) 
1 

1 



±Bl (I? + T* + Tl) + l -B T 2 (TlTl + T 2 T 2 + T 2 T 2 ) 



+ Mj + Y * Tl {T * + T " ] |2 

+ lY^Ti + Y^niTl + TiJf] , 

which has minimum in the VEV direction given by 

Vr 



\Yi T Ti 



Y?T 2 (T? 



T 



(1,1,1). 



(75) 



(76) 



As same as the gauge non-singlet model, we can assume flavor breaking term as perturbation. As the size 
of Vt is at the same order as V, we put 



V T 
Mp~ 



10- 



(77) 



3.2 Quark and Lepton sector 



Due to the Zg symmetry, the effective Yukawa coupling constants accommodate power of e through the 
superpotential 



W c 



( Mr, 



Mr 



(^)o^<W7 + £(t£: 



9( 9i +dj) 



(78) 
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where qi, Ui, di are Zg charge of Qi, Uf, D\ respectively. As the results, the mass matrices of quarks are given 
as follows 



M u ~ e 3 e 2 e ~ 1 1 e diag(e 3 , e 2 , 1) e 1 e \ v' u = L„diag(m u , Clt )i2|„ (79) 
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M d ~ I e 3 e 2 e 2 J ~ I 1 1 e J diag(eV 
from which we get Cabbibo-Kobayashi-Maskawa matrix 

Vckm = L\L d ~ | 1 1 el, 




and quark masses divided by experimental values respectively, 

_jn^_ _ (Al F ) s e*v' u _ rn^_ = (Al F )sfv' u = rn^_ = = 



(m«)exp 1.3 x lO" 3 ' ' (m c ) ex p 0.624 ' ' (m t )exp 173 

m d _ (A^ F )s£ 3 v d = Q ^ m a = {A Y RF ) s e 2 v' d = Q ^ m b = {A Y RF ) s tv' d = J 



(m d ) exp 2.9 x 10- 3 ' (m s ) exp 0.055 (m b )exp 2.89 

(82) 

where e = 0.01, < = 155.3GeV,^ = 77.8GeV, (A^ F ) S = {A B \^. G ) S = 6.647 are used. The renormalization 
factor of top- Yukawa coupling is neglected because it has infrared quasi-fixed point. For the lepton sector, 
the Yukawa coupling constants divided by required values given in Eq.(42) and Eq.(73) are given by 

1 1 1 o nin ^ 



= 1.1, — =2 = — — — =0.19, — =* = ---- =0-16, 



(Y^cxp 0.875 • ' (Y 3 B )exp 5.15 x 10~ 2 ' ' (Y 2 B ) C xp 6.25 x lO" 6 

Y N 1 Y N f Y N 1 

-7.5, -J_=_L- = 0.40 ) _^ = _i_ = 2.7. (83) 



(Y 2 JV )exp 0.134 ■ ' (r 3 Ar )exp 0.025 ' ' (Y 4 N ) e xp 0.370 

Where we used running masses of quarks and charged leptons [12]: 

m u (m z ) = 1.28+p 3g(MeV), m c (m z ) = 624 ± 83(MeV), m t (m z ) = 172.5 ± 3.0(GeV), 

m d (m z ) = 2.91± 1 ;2 o(McV), m s (m z ) = 55±^(MeV), m b (m z ) = 2.89 ± 0.09(GeV), (84) 

m e (m z ) = 0.48657(McV), m^(m z ) = 102.72(MeV), m T (m z ) = 1746(MeV). 

The discrepancies between the estimated values and experimental values in Eq.(82) and Eq.(83) are easily 
recovered by multiplying O(l) coefficients (Y)q. CKM matrix 




(^CifM)exp 0.23 1 4.1 x lO" 2 , [11] (85) 

\ 0.8 x 10~ 2 3.9 x 10~ 2 1 / 

is also recovered from Eq.(81). Therefore our procedure works well in quark and lepton sectors. 
3.3 g-quark sector 

The leading terms of single g-quark interactions are given by 



1 ( X 



Wb = aG[m;I (Yr"MTiGl+T 2 G^+T 3 G c 3 ]U^ b 



q a +l 

M P \ M P J 



+ Tt(tt-) (y® L °°)o[TxGl + TiG\ + T 3 G%]Q a L 3 



M P \Mp) 



(Y a ^ G ) [Vs(T 2 G c 2 - T 3 G§)ii + (T 2 G C 2 + T 3 G C 3 - 2T x G\)L 2 \Q a , (86) 
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where the contribution from E[ D r'° is omitted because p — > t + X is impossible. Note that Y^ G ,Y EUG 
are suppressed by (Vr/Mp) 3 . In the quark mass basis, trilinear coupling matrix and vectors are given as 
follows 



{RlY UDG R d ) ab 
{Y QL ° G L u ) a 
(Y OL « G L u ) a 



1 e e 2 
e 1 e 
e 2 e 1 

lie 
e(e,e,l) | 1 1 e 
eel 



e 4 e 3 e 3 \ 




r 




$ f e : 




1 


;)-( 






1 





e 4 e 3 e 3 
e 3 e 2 e 2 



(87) 



1 1 e 
1 1 e 
e e 1 



^V,e), 
(e,e,l). 



(89) 



As the coupling constants are large enough, the problem of long life time of g-quark is solved. Integrating 
out the scalar g-quarks, would-be the most largest contribution to proton decay is given by 



-p— K° 



MpM G 



($ 3 ) 2 (S c )'(s c )' M '[v / 3ei((T 2 ) 2 - (T 3 ) 2 ) + e 2 ((T 2 ) 2 + (T 3 f - 2 (Ti) 2 )], (90) 



where ei, e 2 are linear combinations of fi and r. Interestingly, this interaction vanishes in the VEV direction 
given in Eq.(76). This means the contributions from three scalar g-quarks are canceled. Therefore the 
dominant contribution to proton decay is given by 



C 



p^re+K° 



6 5 V£ 

M 2 p M% 



from which we get 

T(p -> e+ + A' ) = 

Substituting 

e = 0.01, m N = m p = 940 McV, m B ' 
and Eq.(59), Eq.(64) and Eq.(77), we get 







\ A RF 
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32^/ 2 


_ e5 {m P/ 


1 M%_ 





A RF {u c )'{s c )'u'e, 



m N 
m B > 



(F-D) 



at 



1150 McV, m K o = 498 McV, [11] 



T(p -> e+ + A ) = 1.69 x 10" 



/ lOOOGeV V 
V ) 



GcV. 



For the experimental bound r(p —> e + + K°) > 150 x 10 30 [years] [11], mass bound as follows 

Mq > 1.0 TcV 

must be satisfied. The experimental bound r(p — > v + K + ) > 670 x 10 30 [years] [11] for the operator 



C 



p—u*K+ 



e% 2 
MpM% 



A RF {u c )'{s c )'d'v 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 



gives weaker mass bound as follows 
T{ V ^D + K+) = 



32tt/ 2 ^ 

= 0.20 x 10~ 64 
M G > 0.9 TeV. 



Mp) Ml 
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(97) 
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Note that these bounds should not be taken seriously, because there is 0(10) ambiguity coming from SUSY 
breaking parameter in Vt, this constraint also has such ambiguity. The important point is that we can 
expect the experimental observation of proton decay for TeV scale g-quark in near future. 

Finally we give a short comment about flavor breaking effects on cancellation. Due to the perturbation 
from flavor breaking squared mass terms, scalar g-quark squared mass and squared VEVs of flavons receive 
0{m 2 B / m 2 su SY ) contaminations, where m 2 B is flavor breaking squared mass. As the results, the cancellation 
is spoiled. However, if flavor breaking terms are small enough to satisfy the condition 

TYl 2 

r B = 2 B < e, (99) 

m SUSY 

the cancellation mechanism works effectively. In the opposite case of tb > e, p — > fi + K° dominates the 
proton decay width. Therefore the size of r B affects proton decay channels significantly. There is interesting 
correlation through r B between proton decay channel and degree of degeneracy of scalar g-quark masses. 
Note that too small rs causes the appearance of pseudo-Nambu-Goldstone boson (pNGB). However this 
may be not a serious problem, because if there is the VEV hierarchy such as v s 3> v u ,Vd, which makes Si 
dominant in pNGB, then the interactions of this pNGB with quarks, leptons and weak bosons are very weak. 



4 Dirac neutrino model 

As is shown in previous section, the stability of proton is realized by strong suppression factor of (T) /Alp ~ 
10~ 8 . It is not unnatural to expect that the small neutrino mass is realized by same mechanism. In this 
section, we construct Dirac neutrino model based on 54 flavor symmetry We eliminate $, <f> c and change 
the flavor assignment as given in Table 5. To forbid the bilinear term Hf Li, we change the extra gauge 
symmetry to one linear combination of two extra U(l) gauge symmetries as defined by 

X e = X cos 6 + Z sine, (100) 

and assume RHN and S have non-zero charge of Xg. The value of 9 does not affect our analysis. 
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Table 5: S4 x Z4 x Zg assignment of superfields (Where the index i of the S4 doublets runs i = 1,2, and the 
index a of the S4 triplets runs a = 1, 2, 3.) 

As the quark and charged lepton Yukawa interactions are not modified, we consider only Yukawa inter- 
action of neutrino which is given by 

W N = ^J-(Y 1 JV ) (T 1 7V 1 C + T 2 NZ + TaNSHHYL! + H¥L 2 ) 
X 2 

j^{Y 2 N ) [V3{T 2 N^ - T 3 N^){HYL 2 + H^Lx) + {T 2 N< + T 3 A| - 2Ti.N?){H?Li - H? L 2 )] 
?L(Y 3 N )o[V3(T 2 Ni - T 3 NI)HV + (T 2 N C 2 + T 3 N£ - TT^H^. (101) 
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Substituting the VEVs given in Eq.(74), Eq.(76) and Eq.(77) for X and T a , we get the effective superpotential 
as follows 

W N = Y X N (iff L x + H% L 2 )(NZ + iV 2 c + iV 3 c ) 

+ Y 2 N [V3(NI - NC)(H?L 2 + HU Ll ) + (2JVf - 7V 2 C - N^H?^ - H$ ' L 2 )] 

+ Y 3 N L 3 [V3(iV 3 c - Ni)H? + (27V- - iV 2 c - NDH?} , (102) 

Y&. = (Y&, 3 )oe 2 ( ~ 0(HT 12 ), (103) 



\V3Mpy 

from which Dirac neutrino mass matrix is given by 

((mi + 2m 2 )c„ mic u - m 2 (c u + Visu) m x c u - m 2 (c„ - \/3s u ) 
(m 1 -2m 2 )s u mis u + m 2 (s u - \/3c u ) mis u + m 2 (s„ + V3c u ) ), (104) 
2m 3 s u -m 3 (s u + y/3c u ) m 3 (-s u + V3c u ) 

mi = Y-^Vu, m 2 = Y 2 N v u , m 3 = Y 3 N v u , (105) 

where we can define ui\ and 7713 as real and non-negative and m 2 as complex without loss of generality. In 
the VEV direction 6 U = B c \ = Ob, two large angles of charged lepton and neutrino mixing matrix are canceled 
and make it difficult to realize two large mixing angles of MNS matrix. Therefore we select the condition 
Eq.(17) and put d = and 9 U = j, then the charged lepton mass matrix is given by 

/ mf \ 
Mi = m§ J , (106) 
\ ml j 

from which we get 

/ 1 \ 

V^M l MjV l = diag((ml) 2 ,(ml) 2 7 (m e 1 ) 2 ) = diag(mlml,m 2 T ) 7 = 1 . (107) 

V 1 oj 

To realize experimental results, the conditions given as follows 

V}M* D MlV v = diag{ml^ml 2 ,ml 3 ) = Ml i& g (108) 

= -±=Vi ( 1 1 ) ( 1 I [ -1 V2 



V6 




Vs 



M* D Ml - KM^/J - diag(l, 1, l)m 2 V2 + Am^KdiagC-rv, 0, 1)K + , (HO) 

(111) 



Miag^ =0^1,1,^^71^ 



"'"32 

must be satisfied. Eq.(llO) is rewritten as follows 



(3/2)m 2 + 6|m 2 | 2 -?7i 2 2 (3/2)m 2 6m> 3 

(3/2)m? (3/2)m 2 + 6|m 2 | 2 -m. 2 2 

6m 2 m 3 6to§ — m 2 2 

3-?v 3 + r„ -2?v-3\/2; 

3 + ?v 3-r„ 2r„-3V2A |, (112) 

-2r„ - 3^A* 2r„ - 3\/2A* -4r„ 



':!2 



6 

where 0(A 2 ) terms are neglected. From this equation, we get 



sin0 13 = A = A* = ^r„ = ^^f 1 = 0.014 (0 13 = 0.8°), (113) 



3 " 3 A 



mi = 0.029 eV, m 2 = -0.0033 eV, m 3 = 0.0025 eV 

m Vl = 0.0035 eV, m„ 2 = 0.0094 eV, m„ 3 = 0.051 eV, (114) 
= 2.9 x 10~ 12 , Y 2 N = -0.33 x 10~ 12 , F 3 W = 0.25 x 10~ 12 . (115) 
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The small discrepancies between Eq.(103) and Eq.(115) are recovered by multiplying 0(1) coefficients (Y n )q. 

By the modification of Zg charge of L3, the proton decay width is dominated by p — > e + K° given in 
Eq.(91) because suppression factor is reduced from e 5 to e 4 . From the experimental bound, we get 

1.0 TeV , N 

M G > ^— = 10 TeV. (116) 

v e 

Note that Y NDG is suppressed by (V T /M P ) 2 . 

5 Conclusion 

We have considered proton stability based on £4 symmetric extra U(l) models. Without Froggatt-Nicslen 
mechanism, most stringent bound for proton decay channel is given by r(p — > e + 7r°). As the single quark 
interaction is doubly suppressed by the VEV of gauge non-singlet fiavon, g-quark life time becomes very 
long. Therefore the allowed region for fiavon VEV is very narrow for TeV scale g-quark. 

Introducing Froggatt-Nielsen mechanism, as we can weaken the S4 fiavon VEV suppression, g-quark life 
time becomes short enough. From the naive power counting, we can expect p — > p, + K° would dominate 
the proton decay width, however, corresponding operator vanishes by cancellation and do not contribute to 
proton decay. Therefore our model predicts p — > e + K° dominates the proton decay width. This conclusion 
is not modified in Dirac neutrino model. 
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